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Abstract 

We analyzed the metal distribution of the Cygnus Loop using 14 and 7 point- 
ings observation data obtained by the Suzaku and the XMM-Newton observatories. 
The spectral analysis shows that all the spectra are well fitted by the two-kTg non- 
equilibrium ionization plasma model as shown by the earlier observations. From the 
best-fit parameters of the high-fcTg component, we calculated the emission measures 
about various elements and showed the metal distribution of the ejecta component. 
We found that the distributions of Si and Fe are centered at the southwest of the 
geometric center toward the blow-out region. From the best-fit parameters, we also 
estimated the progenitor mass of the Cygnus Loop from our field of view and the 
metal rich region with a radius of 25 arcmin from the metal center. The result from 
the metal circle is similar to that from our entire FOV, which suggests the mixing of 
the metal. From the results, we estimated the mass of the progenitor star at 12-15M0. 

Key words: ISM: abundances — ISM: individual (Cygnus Loop) — supernova 
remnants — X-rays: ISM 

1. Introduction 

The supernova (SN) explosions blow out the various heavy elements generated by the 
nucleosynthesis process inside the progenitor stars. Meanwhile, the blast wave formed by the 
SN explosion sweeps up and heats the interstellar matter (ISM). In this way, the abundance of 
the progenitor star or the ISM are measured from the observation of the supernova remnant 
(SNR). 

The Cygnus Loop is one of the brightest SNRs in the X-ray sky. Its age is estimated 
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to be ~10,000 yr and the distance is comparatively close to us (540 pc; Blair et al. 2005). The 
large apparent size (2°. 5 x 3°. 5; Levenson et al. 1997) enables us to study the plasma structure 
of the Loop. 

Although the Cygnus Loop is an evolved SNR, a hot plasma is still confined inside 
the Loop (Tsunemi et al. 1988; Hatsukade & Tsunemi 1990). Miyata et al. (1998) observed 
the Loop with the Advanced Satellite for Cosmology and Astrophysics (ASCA), and detected 
strong highly-ionized Si-K, S-K, and Fe-L lines near the center of the Cygnus Loop. They 
concluded that a hot plasma, a "fossil" of the supernova explosion, is left in the core of the Loop. 
Tsunemi et al. (2007) observed the Cygnus Loop along the diameter from the northeast (NE) 
to the southwest (SW) with XMM-Newton and studied the radial plasma structure. From the 
spectral analysis, they showed that the Cygnus Loop consists of two components with different 
temperatures and metal abundances. They concluded that the low-ZcTg component originating 
from the cavity-wall component surrounds the high-ZcTg ejecta component. In addition, they 
measured the metal abundances of the high-kT^. component and showed the metal distribution 
of the ejecta. The results indicate that the abundances are relatively high (~5 times solar) and 
each element is non-uniformly distributed: Si, S and Fe are concentrated in the inner region 
while the other elements such as O, Ne and Mg are abundant in the outer region. Katsuda 
et al. (2008a) and Kimura et al. (2009)^ expanded the observations southward and northward 
respectively with the Suzaku observatory in 7 and 10 pointings, and examined the plasma 
structure in their Field of View (FOV). Katsuda et al. (2008a) divided their FOV into 119 cells 
and, performed the spectral analysis. They showed that Si and Fe are more concentrated in 
the south part than that in the north part. This indicates a clear asymmetric structure of the 
metal abundances. Kimura et al. (2009) divided their FOV into 45 rectangular regions from 
NE to SW and also showed the asymmetric distribution of each heavy element; the ejecta of 
O, Ne, Mg are distributed more in the NE, while those of Si, Fe are distributed more in the 
SW. Tsunemi et al. (2007) and Kimura et al. (2009) also calculated the progenitor mass of the 
Cygnus Loop and concluded that the Cygnus Loop is originated from the 12-15MQexplosion. 
Levenson et al. (1998) also estimated the progenitor mass to be 15M0from the size of the cavity. 
These results show that the progenitor star of the Cygnus Loop was a massive star and caused 
a core-collapse explosion. It is striking that an asymmetric explosion is suggested as the origin 
of the Cygnus Loop from the metal distributions. 

The Cygnus Loop is a typical shell-like SNR; this structure is thought to be generated 
by a cavity explosion (Levenson et al. 1997). The morphology is almost circular while we can 
see the breakout in the south called the "blow-out" region (Aschenbach & Leahy 1999). The 
origin of the "blow-out" is not well understood. Aschenbach & Leahy (1999) have explained 
this extended structure as a breakout into a lower density ISM. On the other hand, Uyaniker 
et al. (2002) suggested the existence of a secondary SNR (named G72. 9-9.0) in the south based 

1 The data available at http://arxiv.org/pdf/0810.4704vl 
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on radio observation and some other radio observations support this conclusion (Uyaniker et al. 
2004; Sun et al. 2006). Recently, Uchida et al. (2008) observed the blow-out region with the 
XMM-Newton observatory and showed that the X-ray spectra of this region consist of a two-kTe 
plasma component. Judging from the plasma temperatures and the metal distributions, Uchida 
et al. (2008) concluded that the X-ray emission is consistent with a Cygnus Loop origin and 
that the high- and \ow-kTe components derived from the ejecta and the cavity-wall, respectively. 
They also showed that the X-ray shell is thin in the blow-out region and concluded that the 
origin of the blow-out can be explained as a breakout into a lower density ISM as proposed by 
Aschenbach & Leahy (1999). 

In this paper, we used 14 and 7 pointings observation data obtained by the Suzaku and 
the XMM-Newton observatories. We reanalyzed all the data to study the metal distributions 
inside the Cygnus Loop. 

2. Observations 

We summarized the 21 observations in table 1. All of the Suzaku data were analyzed 
with version 6.5 of the HEAsoft tools. For the reduction of the Suzaku data, we used the version 
9 of the Suzaku Software. The calibration database (CALDB) used was the one updated in July 
2008. We used the revision 2.2 of the cleaned event data and combined the 3x3 and 5x5 event 
files. The P03, P04, P05, P06, P07, P09, Pll, and P20 data were taken by using the spaced 
row charge injection (SCI) method (Prigozhin et al. 2008) which reduces the effect of radiation 
damage of the XIS and recovers the energy resolution. In order to exclude the background flare 
events, we obtained the good time intervals (GTIs) by including only times at which the count 
rates are within ±2(j of the mean count rates. 

Since the Cygnus Loop is a large diffuse source and our FOV are filled with the SNR's 
emission, we cannot obtain the background spectra from our FOV. We also have no background 
data from the neighborhood of the Cygnus Loop. We therefore applied the Lockman Hole data 
for the background subtraction. We reviewed the effect of the galactic ridge X-ray emission 
(GRXE). The flux of the GRXE at / = 62°, \h\ < 0°.4 is 6 x lO^^^ g^g/cniVs/deg^ (0.7-2.0 keV) 
(Sugizaki et al. 2001). Although the Cygnus Loop (/ = 74°, h = —8°. 6) is located outside 
of the FOV of Sugizaki et al. (2001), this value gives us an upper limit of the GRXE at 
the Cygnus Loop. Meanwhile, the total count rate of the Cygnus Loop is 1659 counts/s (0.8- 
2.01 keV) (Aschenbach & Leahy 1999), the flux is estimated to be 1.68 x 10~^° erg/cm^/s/deg^, 
assuming that the effective area of the ROSAT HRI is 80 cm^ and that the Cygnus Loop is a 
circle 3°.0 in diameter. This value is consistent with the results from our FOV. Therefore, we 
concluded that the effect of the GRXE on the Cygnus Loop is vanishingly small. The solar 
wind charge exchange (SWCX) is also considered to be one of the correlates of the soft X-ray 
background below IkeV (Fujimoto et al. 2007). However, in terms of the Cygnus Loop, we 
consider that the SWCX is negligible because of the prominent surface brightness of the Loop. 
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Thus, the Lockman Hole data obtained in 2006 and 2007 were apphed for the background 
subtraction. We selected data whose observation dates were close to those of the Cygnus Loop 
observations. Since there were no photons above 3.0 keV after the background subtraction, the 
energy ranges of 0.3-3.0 keV and 0.4-3.0 keV were used for XISl (back-illuminated CCD; BI 
CCD) and XIS0,2,3 (front-illuminated CCD; FI CCD), respectively (Koyama et al. 2007). 

All of the XMM-Newton data were processed with version 7.1.0 of the XMM Science 
Analysis System (SAS). The current calibration files (CCFs) used were the one updated on 
June 2008. We used the data obtained with the EPIC MOS and pn cameras. These data were 
taken by using the medium filters and the prime full-window mode. We selected X-ray events 
corresponding to patterns 0-12 and flag = for MOS 1 and 2, patterns 0-4 and flag = for pn, 
respectively. In order to exclude the background flare events, we determined the GTIs in the 
same way as the Suzaku data. After filtering the data, they were vignetting-corrected using the 
SAS task evigweight. For the background subtraction, we employed a blank-sky observations 
prepared by Read & Ponman (2003) for similar reason with the Suzaku^s case. After the 
background subtraction, the energy ranges of 0.3-3.0 keV were used for each instrument. 

3. Spectral Analysis 

To investigate the plasma structure of the Cygnus Loop, we divided the entire FOV 
into several box regions outlined in white line in Fig. 1. In order to equalize the statistics, 
we initially divided all images of XISl or M0S2 into two parts and if each divided region has 
more than 10,000 photons, it was once again divided. In this way, we obtained 415 box regions. 
Each region contains 5,000-10,000 photons for XISl and M0S2. The side length of each box 
ranges from 2.2 to 14arcmin. Therefore box sizes are not smaller than the angular resolution 
capability of the Suzaku XIS. We grouped 415 spectra into bins with a minimum of 20 counts 
so that statistics are appropriate. From the earlier observation of the NE to the SW regions 
along the diameter, Tsunemi et al. (2007) showed the plasma structure of the Cygnus Loop as 
follows: the high-/cTe ejecta component is surrounded by the low-kT^ ISM component. They 
found that the spectra from the inner regions of the Cygnus Loop consist of the two-component 
non-equilibrium ionization plasma. In order to examine the metal distributions of the ejecta of 
the Cygnus Loop, we need to separate the X-ray spectra into the ejecta and the surrounding 
cavity-material component. Then, we fitted 415 spectra extracted from our FOV by the two-kTe 
non-equilibrium ionization (NEI) plasma model. We employed wabs (Morrison & McCammon 
1983) and vnei (NEI ver.2.0; Borkowski et al. 2001) in XSPEC version 12.4.0 (Arnaud 1996). 
In order to generate a response matrix file (RMF) and an ancillary response file (ARE), we 
employed xisrmfgen (Ishisaki et al. 2007) and xissimarfgen for the Suzaku data, rmfgen and 
arfgen for the XMM-Newton data. We generated ARFs of the Suzaku data for "extended 
sources" . 

First, we set all the abundances of the heavy elements free, and found that this model 
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Table 1. Summary of the 21 observations 



Obs. ID 


Obs. Date 


Coordinate (RA, DEC) 


Position Angle 


Effective Exposure 


Suzaku Observations 


501014010 (P3) 


2007-11-14 


20''52'"09 .9, 


31 36 43 


4 


58°. 4 


7.5 fcsec 


501015010 (P4) 


2007-11-14 


20^51 11 .8, 


31 22 08 


4 


58°. 6 


18.3 ksec 


501016010 (P5) 


2007-11-15 


20"50 11 .3, 


01 ni A r\ f 1 

31 10 48 





58°. 8 


19.3 ksec 


501017010 (P6) 


2007-11-11 


20''49 11 .3, 


30 59 27 


6 


59°. 


28.7 ksec 


501018010 (P7) 


2007-11-12 


20"48 18 .7, 


30 46 33 


6 


59°. 2 


21.0 ksec 


501019010 (P9) 


2007-11-12 


20"47 14^.2, 


30 36 10 


8 


59°. 4 


16.2 ksec 


0UoUi30UiU (,i llj 


onnc nc; no 
zUUo-Uo-U9 


zU 4y 4o 


61 oU io 


U 


00 .0 


zz.z Ksec 


501029010 (P12) 


2006-05-09 


20''55"00".0, 


31°15'46" 


8 


58°. 7 


13.2 ksec 


501030010 (P13) 


2006-05-10 


20''53"59".3, 


31°03'39" 


6 


58°. 9 


13.9 ksec 


501031010 (P14) 


2006-05-12 


20''52"58".8, 


30°51'32" 


4 


59°. 1 


18.2 ksec 


501032010 (P15) 


2006-05-25 


20''51"58".6, 


30°39'10" 


8 


59°. 3 


17.4 ksec 


501033010 (P16) 


2006-05-22 


20''50"58".8, 


30°27'00" 





59°. 5 


20.0 ksec 


501034010 (P17) 


2006-05-22 


20''48"49".7, 


30°00'21" 


6 


60°. 


13.9 ksec 


503056010 (P20) 


2008-05-10 


20''48"00".0, 


31°10'30" 





58°. 8 


22.5 ksec 


XMM-Newton Observations 


0082540201 (Pos-2) 


2002-12-03 


20''54'"07".2, 


31°30'51" 


1 


241°. 7 


14.4 ksec 


0082540301 (Pos-3) 


2002-12-05 


20''52"51M, 


31°15'25" 


7 


241°. 7 


11.6 ksec 


0082540401 (Pos-4) 


2002-12-07 


20''51"34".7, 


31°00'00" 





241°. 7 


4.9 ksec 


0082540501 (Pos-5) 


2002-12-09 


20''50"18".4, 


30°44'34" 


3 


231°.4 


12.6 ksec 


0082540601 (Pos-6) 


2002-12-11 


20''49"02".0, 


30°29'08" 


6 


241°. 7 


11.5 ksec 


0405490101 (Pos-8) 


2006-05-13 


20''50"32".2, 


30°11'00" 





69°. 9 


6.5 ksec 


0405490201 (Pos-9) 


2006-05-13 


20''49"54".2, 


29°42'25" 





69°. 8 


3.6 ksec 



5 



cannot reach the physically meaningful results. Therefore, in the low-ZcTe component, we fixed 
the metal abundances to the values from the result of Tsunemi et al. (2007). The abundances of 
the low-kTe ISM component shown by them are as follows: C=0.27, N=0.10, 0=0.11, Ne=0.21, 
Mg=0.17, Si=0.34, S=0.17, Fe(=Ni)=0.20. In addition, they fixed the other elements to the 
solar values. Then, we used these values as the abundances of the emission from the ISM origin. 
As for the high-fcTe component, the abundances of O, Ne, Mg, Si, and Fe were free while we 
set the abundances of C and N equal to O, S equal to Si, Ni equal to Fe, and other elements 
fixed to their solar values (Anders & Grevesse 1989). Other parameters were all free such as 
the electron temperature kT^., the ionization timescale r (a product of the electron density and 
the elapsed time after the shock heating), and the emission measure (EM = / neUndl, where Ue 
and riH are the number densities of hydrogen and electrons and dl is the plasma depth). We 
also set the column density A'h free. 

The spectra are reasonably well fitted by the two-ZcTg VNEI model for almost all regions: 
the values of the reduced shows around ~1.5 and the degrees of freedom are 300-400. Some 
regions overlap each other and we fitted these spectra separately. In cases like this, we compared 
each best- fit parameter values and confirmed that they are just within the uncertainties. The 
example spectra and the best-fit curves are shown in figure 2 and the best-fit parameters are 
shown in table 2. These two spectra are taken from the region where Fe and Si are the most 
abundant or the most depleted, respectively. The extracted regions are shown in figure 1. 

Figure 3 shows the electron temperature distribution of each component. The left panel 
shows the best-fit kTe parameters of the high-/cTe component that represents the plasma tem- 
perature of the ejecta component. From the left panel, we show the inhomogeneous temperature 
distribution of the ejecta component. The values of kT^ range from 0.4keV to 0.7keV and they 
are lower in the SW part than those in the NE part. The ejecta temperature of the Cygnus 
Loop shows smooth decrease from NE to SW. The averaged temperature is ~0.49keV. The 
right panel shows the map of the kTe parameter obtained from the \ow-kTe component that 
represents the temperature of the swept-up cavity material. The averaged value is ~0.23keV 
and it ranges from 0.2 keV to 0.3 keV. The temperature of the low-kT^ component near the 
center of the Loop is lower than that of the surrounding temperature. From figure 3, we clearly 
separated the high-kTg component and the low-kT^. component just as the observation obtained 
in Tsunemi et al. (2007), Katsuda et al. (2008a), and Kimura et al. (2009). 

Figure 4 shows the EM distributions of the heavy elements. The images are smoothed by 
Gaussian kernel of o"=2.8 arcmin. The color code scale is normalized by the maximum values. 
From figure 4, the EM of Mg is relatively low compared with those of the other elements. The 
EM values of O, Ne, Mg are higher in the NE than those around the center while the EM values 
of Si and Fe are larger in the center and decrease outward. These results are consistent with 
those of Tsunemi et al. (2007), Katsuda et al. (2008a), and Kimura et al. (2009). These earlier 
studies were all conducted from NE to SW in their FOV. Our analysis revealed the plasma 
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Fig. 1. Left panel shows the ROSAT HRI image of the entire Cygnus Loop overlaid with the spectral 
extraction regions with white rectangles. Right two panels show the extraction regions for Suzaku (top 
panel) and XMM-Newton (bottom panel). 

structure from the northwest (NW) to the southeast (SE) of the Loop. The averaged value 
of each EM [xlO^^ cm-^] is 3.20 (0=C=N), 0.39 (Ne), 0.08 (Mg), 0.81 (Si=S), 0.59 (Fe=Ni), 
respectively. We show the geometric center estimated by Levenson et al. (1998) with the black 
X-mark in figure 4. It is located on a = 20^51^2^,5 = 31°0l'37" (J2000) which is determined 
by fitting the ROSAT HRI map of the Loop by the model circle. 

4. Discussion 

4-1. Low abundances of the low-kTf. component 

We analyzed all the spectra with using a two-component model: the low-kTe compo- 
nent forms a shell-like structure and the high-fcTg component fills inside (Tsunemi et al. 2007; 
Katsuda et al. 2008a; Kimura et al. 2009). The fixed abundances of the low-kTe component are 
depleted compared to the solar values. Cartledge et al. (2004) measured the abundance of the 
interstellar oxygen in the solar neighborhood along 36 sight lines and showed 0.4 times the solar 
value. Wilms et al. (2000) employed 0.6 of the total interstellar abundances for the ISM oxygen 
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Energy (keV) ' Energy (keV) 



Fig. 2. Example spectra from the regions where Fe and Si are the most abundant (region A: left) and 
the most depleted (region B: right), respectively. The best-fit curves are shown with solid line and two 
components are shown with dotted lines. The high- and the low-fcTe components cross at ~0.65keV 
(region A) and '^0.7 kcV (region B), respectively. The residuals arc shown in the lower panels. In the 
left panel, black, red, green, blue correspond to the XIS 0, 1, 2, 3. In the right panel, black, red, green 
correspond to the XIS 0, 1, 3. 




0.1 0.2 0.3 0.4 0.5 0,6 0.7 



Fig. 3. Electron temperature distribution of each component. The left and right panel shows the distri- 
bution of the high- and low-fcTe component, respectively. The values of fcTg are in units of kcV. 
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Fe 



1 2 

Fig. 4. EM distributions of the heavy elements in ?hc logarithmic scales. The black X-mark shows the 
geometric center of the Loop (Levenson et al. 1998). The blue contour shows the EM level of 1/e of the 
maximum around the geometric center. The light blue cross-mark and the dotted circle represent the 
"metal center" and the "metal circle", respectively (see 4.2). The values are in units of lO^^cm"^. 



Table 2. Spectral fit parameters 



Parameter 


region A 


region B 


iNfj [iU cm J 


6.0 ± U.z 


Z.D it O.Z 


Low-fcre component: 








n 1 n _i_ n m 
u.iy ± U.Ul 


O.zo ± U.Oi 




U.z < 


^nxecl j 


AT 
IN 


U.iU 


(tixed j 


U 


nil 
U.ii 


(tixed j 


AT^ 

IN 6 


n o 1 


(nxedj 


iVig 




(^nxecL j 


bi 


0.64: 


(nxedj 


o 
t) 


n 1 7 


(nxed ) 


^ AT'N 

-be(=iNij 


A on 


(nxedj 


log r 


n /in +0.05 


1119 +0.03 
rr-rz _o.o3 


T7\ i\/r finis — F^^ 

hM [10 cm °\ 


o o 1 n 1 n 

6.23 ± 0.19 


A 0^ 1 n 1 H 

4.27 ± 0.14 


High-fcTe component: 






fcie [keVJ 


r\ /t o 1 Am 

(J.4z ± U.Ul 


A 1 A A1 

0.53 ± (J.Ul 


/ AT^ 


A /I T _i_ n no 
U.4/ ± U.Uo 


A 1 _l_ A AO 

[).6i ± U.Uz 


AT^ 


A Tc: 1 A nn 

0.75 ± 0.09 


A /I T 1 A A C: 

0.47 ± 0.05 


Mg 


0.21 ± 0.10 


0.41 ± 0.08 


bi[~b ) 


O.OD ± 0.4» 


0.49 it O.lo 


Fe(==:Ni) 


2.92 ± 0.05 


0.55 ± 0.03 


log r 


11 7f; _L +0.10 


10.82 ± 0.02 


EM [IQi^cm-S] 


0.55 ± 0.01 


0.68 ± 0.02 


xVdof 


441/308 


527/402 



abundance. The abundances we fixed are significantly lower than those results and it needs to 
be discussed. The metal deficiency of the ISM component of the Cygnus Loop has been shown 
by the earlier X-ray observations (Miyata & Tsunemi 1999). The other observations in X-ray 
by using the Suzaku, XMM-Newton, and Chandra also showed the metal depletion at the rim 
(Leahy 2004; Tsunemi et al. 2007; Katsuda et al. 2008b). Meanwhile, Katsuda et al. (2008b) 
found an enhanced abundance region in the very edge of the NE rim. The abundance there is 
about half the solar value which is much higher than those other rim regions. This result is 
consistent with the observations of the ISM abundance in the solar neighborhood (Wilms et al. 
2000; Cartledge et al. 2004). Therefore, the metal abundances of the periphery of the Loop may 
be higher than those observed in X-ray. Any line of sight through the remnant almost certainly 
intersects various emitting regions. However, the line of sight near the limb can intersect the 
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very limited plasma condition. Therefore, we can expect that the very rim region Katsuda 
et al. (2008b) observed consists of thermal plasma with ISM abundance while other regions 
may include either thin thermal plasma with low abundance or non-thermal emission. Katsuda 
et al. (2008c) tried to fit the spectrum with a combination of thin thermal plasma with ISM 
abundance and the non-thermal emission. They reached no clear result probably due to the 
simplicity of the model. Although most of the rim regions show low abundance from the X-ray 
data analysis point of view, it still remains an open question. 

4.2. Displacement of the metal distribution from the geometric center 

We noticed that the strong EM regions concentrate just south of the geometric center. 
Si and Fe concentrate in the inner part while Ne and O have strong regions both in the inner 
part and the outer part within our FOV. Mg has strong regions only in the outer part. In this 
paper, we will focus on the strong regions near the geometric center. 

The contours in figure 4 show the EM level of 1/e of the maximum near the geometric 
center. The concentration of the inner part of various metals is obviously away from the 
geometric center. The spread of the strong regions indicated by the contour is not clear in the SE 
region due to the lack of observation. Therefore, it is difficult to precisely determine the spread. 
Taking into account the lack of information in the SE region, we speculate that the center of 
the metal distribution is determined by the centroid of the Fe distribution at the south of the 
geometric center. We call it the "metal center" , which is located on a = 
(J2000), as shown by the light blue cross- mark in figure 4. 

We found that the metal center is separated from the geometric center by 25 arcmin 
toward the south. It is not clear why it is displaced from the geometric center. On the simple 
assumption that the SN explosion occurred symmetrically in the uniform ISM, the metals should 
distribute symmetrically and its center should overlap with the geometric center. The Cygnus 
Loop shows a nearly circular shape with an exception of the south blow-out region. However, 
the X-ray surface brightness and the optical filamentation suggest that the ambient medium of 
the Cygnus Loop is not uniform. We consider that its non-uniformity is prominent toward the 
south blow-out region. If the south blow-out is originated from the lower-density cavity wall or 
the ISM, the displacement of the metal distributions to the south could be naturally explained. 
In general, the non- uniformity produces the asymmetrical reverse shocks. Taking into account 
the estimated age of the Loop, the reverse shocks should have reached the center of the Loop. 
However, the times for the reverse shocks to return to the center, significantly depend on 
the ambient density in the direction which the forward shock propagates. From the recent 
theoretical calculation, Ferreira & de Jager (2008) obtained the relation of oc (pism)"^^'^, 
where pisM represents the ambient density. Therefore, in directions where pisM is higher, the 
forward shock decelerates and the reverse shock propagates toward the center more quickly. 
From the morphological point of view, it is suggested that the cavity wall must be very thin in 
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the south blow-out region, ahhough the origin of the blow-out is still not clear. 

As described in section 1, from the radio observations, Uyaniker et al. (2002; 2004) and 
Sun et al. (2006) claim that the extra SNR exists in the south and interacts with the Cygnus 
Loop. Their main arguments are the difference of the radio morphology and the polarization 
intensity between the main part of the Loop and the south blow-out. However, based on the 
X-ray data, Uchida et al. (2008) obtained no evidence that a smaller SNR exists at the same 
distance to the Cygnus Loop. They showed that the X-ray spectra in the south are interpretable 
by the Cygnus Loop origin. They also found that the emission from the cavity material of the 
Loop is relatively weak. It suggests that the X-ray shell is very thin in the south. On the other 
other hand, the ambient medium in the north part of the Cygnus Loop is thought to be denser 
than that in the south. The north areas have the strong emission in X-ray and it suggests that 
the forward shock has decelerated enough to become radiative in places. Then, the reverse 
shocks in some part other than south should be formed earlier and reaches the center faster 
than that in the south. In that case, the pressure exerted by the north reverse shock is higher 
than that exerted by the south one. If the metal distributions are formed by the imbalance of 
the pressure surrounding them, they should be created southward from the geometric center. 
Our result is consistent with the assumption that the south blow-out is originated not from an 
extra SNR but from a non-uniformity of the cavity wall. 

4-3. Metal circle and the progenitor star mass 

Although the metal abundance is not uniform within the high-fcTg component, the con- 
tents in our FOV reflect the abundance of the progenitor star. We calculate the metal abun- 
dances around the metal center and in our FOV relative to O in order to compare them with 
the theoretical models. We defined the "metal circle" which is centered in the metal center 
with radius of 25 arcmin. We added the EM of ejecta component of each region weighting each 
area of the region for the heavy elements, O, Ne, Mg, Si, and Fe. These values reflect the total 
amount of each heavy element in our FOV. Then, we obtained the number ratios of each element 
relative to O. Figure 5 shows the number ratios of Ne, Mg, Si and Fe relative to O of the ejecta 
component. The vertical bars show the uncertainties of the values. The black lines represent 
the results from the metal circle and the entire FOV. We also plotted the core-collapse models 
(Woosley & Weaver 1995) for various progenitor masses and Type la SN models (Iwamoto 
et al. 1999) for comparison. The result from the metal circle is similar to that from our entire 
FOV, which suggests the mixing of the metal. We found that the calculated Ne/0, Mg/0, and 
Si/0 are in good agreement with the core-collapse model with the progenitor mass of 12Mq. 
However, the Fe/0 ratio is 10 times higher than that of the 12MQmodel, and it does not fit 
any other models. This may be partly because our FOV does not cover the entire region of 
the Loop. From figure 4, we found that the EM values of O, Ne, Mg are higher around our 
FOV rather than those near the center of the Loop. By contrast. Si and Fe decrease toward 
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Fig. 5. Number ratios of Ne, Mg, Si, and Fe relative to O of the high-fcTg component estimated from 
the metal circle and that from the entire FOV (black lines). Dotted red lines represent the CDDl and 
W7 Type la supernova models of Iwamoto et al. (1999). Dotted blue, light blue, magenta, and green 
lines represent core-collapse models with progenitor masses of 12, 13, 15, 20Mq, respectively (Woosley & 
Weaver 1995). 

the outside of the Loop. If these trend continue to the outside of our FOV, we expect that 
the total EMs of O, Ne, Mg will increase considerably more than those of Si and Fe. In that 
case, Si/0 and Fe/0 estimated by the whole remnant will decrease and approach the values 
the core-collapse models rather than those of the Type la SN models. We can conclude that 
the progenitor star mass of the Cygnus Loop is 12-15Mq. 

5. Conclusion 



We analyzed the metal distribution of the Cygnus Loop using 14 and 7 pointings obser- 
vation data obtained by the Suzaku and the XMM-Newton observatories. We fit all the spectra 
by the two-ZcTg non-equilibrium ionization plasma model (VNEI model) as shown by the earlier 
observations (Tsunemi et al. 2007; Katsuda et al. 2008a; Kimura et al. 2009). The results 
indicate that Si and Fe are more abundant near the center of the Loop than those around our 
FOV. We measured a center of the Si and Fe distributions from the metal spread and called it 
as "metal center" . We found that the metal center is located at the southwest of the geometric 
center toward the blow-out region. 

From the best-fit parameters, we also estimated the progenitor mass of the Cygnus Loop. 
We calculated the number ratios of the heavy elements (relative to O) both from the entire 
FOV and inside the metal circle. The result from the metal circle is not different from that 
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of our entire FOV, which suggests mixing of the metals. The results show that Ne/0, Mg/0, 
and Si/0 are well fitted by the core-collapse models with progenitor masses of 12, 13, 15Mq. 
However, Fe/0 is 10 times higher than those of the models. This may be partly because our 
FOV does not cover the entire region of the Loop. We can conclude that the progenitor mass 
of the Cygnus Loop is 12-15M0. 
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